Lack of Evidence for a Singlet Crystal Field Ground State in the Tb2Ti207 Magnetic Pyrochlore 
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We present new high resolution inelastic neutron scattering data on the candidate spin liquid Tb2Ti207. We 
find that there is no evidence for a zero field splitting of the ground state doublet within the 0.2 K resolution 
of the instrument. This result contrasts with a pair of recent works on Tb2Ti207 claiming that the spin liquid 
behavior can be attributed to a 2 K split singlet- singlet single-ion spectrum at low energies. We also reconsider 
the entropy argument presented in Chapuis et al. as further evidence of a singlet-singlet crystal field spectrum. 
We arrive at the conclusion that estimates of the low temperature residual entropy drawn from heat capacity 
measurements are a poor guide to the single ion spectrum without understanding the nature of the correlations. 

PACS numbers: 75.10.Dg 75.10.Jm 75.40.Cx 75.40.Gb 



Introduction — In some magnetic systems, the lattice ge- 
ometry or the competition between different interactions can 
dramatically inhibit, or frustrate, the development of long 
range order. The failure of a frustrated magnetic system to 
exhibit magnetic order down to zero temperature, giving rise 
to a so-called spin liquid state, is one of the most sought af- 
ter phenomena among strongly interacting condensed matter 
systems. Despite two decades of experimental searches, the 
number of candidate materials that display a spin liquid state 
remains small The Tb2Ti207 insulating material, where 
magnetic Tb"^"*" ions sit on a pyrochlore lattice of corner- 
sharing tetrahedra, is one of these candidates fl. Despite a 
Curie- Weiss temperature, 6'cw ~ — 14 K set by the magnetic 
interactions Jst], Tb2Ti207 does not develop long range order 
down to at least 50 mK The microscopic mechanism by 

which Tb2Ti207 fails to develop long range order at a temper- 
ature scale of approximately 1 K IHIslS] is not understood 10] ■ 
Compounding the difficulty in understanding why Tb2Ti207 
does not order, one notes that Tb2Sn207, seemingly closely 
related at the microscopic level to Tb2Ti207, develops long 
range order at 0.87 K 1113]. 

Since Tb'^^ is an even electron system (electronic config- 
uration ''Fg, L = 3, 5* = 3, J = 6), the existence of 
a magnetic ground state for an an isolated (i.e. assumed 
non-interacting) Tb^+ ion in Tb2Ti207 is not guaranteed by 
Kramers' theorem jjst]. To investigate whether Tb^+ is mag- 
netic, the so-called single-ion crystal field (CF) problem must 
first be solved jUj]. For example, for a perfect cubic ionic 
environment, theory predicts that Tb'^^ would either have a 
sing let or non-magnetic doublet single-ion CF ground state 
llioll . However, in Tb2Ti207, with its Fd3m symmetry, the 
Tb'^"'" environment displays a very large trigonal distortion 
away from cubic symmetry ]|3|. Early investigations found 
that this distortion endows Tb'^+ with a magnetic CF doublet 
characterized by two mutually time-reversed conjugate states. 



and I -00") iIjI. The states, |V'o^), are such that all ma- 
trix elements of the raising and lowering angular momentum 
operator, ,J^, vanish while (V'^|J^|V'^) = ±|(J^)| ~ 3.4 
m- Since (V'^I^^IV'^) = ±|(J'')|^m,^' *e Tb^^+ moment 
within its CF ground state can be described by a classical Ising 
spin with a moment that points "in" or "out" of the reference 
primitive tetrahedral unit cell to which it belongs. This makes 
Tb2Ti207 a relative of the Ising spin ice compounds lHQ]- 

The lowest excited CF state is also a doublet, I^A^ ), at an 
energy approximately 1.6 meV ^ 18 K above the ground dou- 
blet 13, 1 l.jj. Recent theoretical work has argued that the 
proximity of Tb2Ti207 to the zero temperature transition from 
an "all-in/all-out" q = Neel to a spin ice state llsf], along with 
the exchange and dipole-dipole interaction-induced admixing 
of l?/!^) and l"!/"^), constitute two key ingredients as to why 
Tb2Ti207 does not develop long range order. 

In conventional (unfrustrated) magnets, spin-lattice cou- 
plings typically play an insignificant role in the develop- 
ment of long range order. In highly frustrated magnetic sys- 
tems, however, spin-lattice couplings can lead to a combined 
magnetic-lattice ("spin-Peierls") transition to long range mag- 
netic order that reduces the magnetic frustration [Il3l . as 
observed in the highly frustrated antiferromagnet ZnCr204 
spinel compound lll4ll . In Tb2Ti207, Mamsurova and co- 
workers long ago reported an unusually large anomalous field - 
dependent thermal expansion, indicating an important spin- 
lattice coupling in this material lIlSll . In more recent x-ray 
diffraction experiments. Ruff et al. found evidence for a ten- 
dency of the Tb2Ti207 lattice to undergo a cubic tetragonal 
deformation but, down to 300 mK and in zero field, no equi- 
librium cubic to tetragonal transition was observed iflHI . Only 
in the presence of very high magnetic fields (^ 29 T), does the 
system show any evidence for such a structural phase transi- 
tion ITtIi . In contrast, Chapuis and co-workers very recently 
argued that a tetragonal deformation does actually exist in 
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Tb2Ti207 in zero field that splits the above I'iI'q) doublet into 
two non-magnetic singlets separated by an energy scale, S, 
with S « 1.8 K IT^ . They invoke previously published in- 
elastic neutron scattering (INS) data 1I12I1 to suggest that such a 
large splitting of the ground doublet compared to the exchange 
and dipolar interactions between Tb^^^ ions is responsible for 
inhibiting the spontaneous development of long range order. 
Reference llSlI also present data for the temperature depen- 
dence of the magnetic entropy, S{T), which, they claim, by 
falling below Rln 2 at low temperature, further supports the 
evidence for a split doublet. Very recently, Bonville and col- 
laborators lll9ll have built further on Chapuis et al. 's split dou- 
blet picture to advocate that the failure of Tb2Ti207 to order 
is due to the sub-critical value of the interactions compared to 
the singlet-singlet gap 5. 

In this paper, we argue that the evidence for a split doublet 
of energy scale as large as (5 ~ 1.8 K in Tb2Ti207 as proposed 
m Refs. iS [il is not compelling. Using new high resolu- 
tion INS data, we show that there is no evidence for a split 
doublet in this material with an energy splitting greater than 



0.2 K, a factor 10 or so smaller than the proposed fllSL 11911 
singlet-singlet gap 5. Secondly, we argue that by neglecting 
correlations that develop in the collective paramagnetic (spin 
liquid) phase of Tb2Ti207, the authors of Ref. ifisll are in prin- 
ciple unable to draw any conclusion about the nature of the CF 
state of Tb'^+ in Tb2Ti207 on the basis of a measurement of 
S{T). We illustrate that point via the calculation of S{T) for 
a toy model which, while possessing a ground state doublet 
and lacking a transition to long range order, does display an 
S{T) that falls below Rln 2 at low temperature. 

Inelastic neutron scattering results — We first proceed to 
show that a singlet-singlet gap 5 ^ 1.8 K is inconsistent with 
high energy-resolution inelastic neutron scattering data and 
that, in fact, Tb2Ti207 displays quasielastic magnetic spectral 
weight down to energies of at least 0.02 meV, approximately 
an order of magnitude lower in energy than the value (5 ~ 1.8 
K reported in Refs. HHH. 

Time-of-flight neutron scattering data was obtained on sin- 
gle crystal Tb2Ti207 in zero and finite magnetic field applied 
along the [110] (vertical) direction at T=0.4 K. These measure- 
ments used the Disk Chopper Spectrometer (DCS) at NIST 
and employed incident neutrons of wavelength A=4.8 A re- 
sulting in an energy resolution of 0. 1 me V. The H = 0, zero 
magnetic field data set, after integration along the [HHO] di- 
rection, perpendicular to the (OOL) direction and within the 
horizontal scattering plane, is plotted in Fig. 1 a). The quasi- 
elastic magnetic scattering peaks at Q=002 (Q=1.2 A"^). Pre- 
vious work (3, Eltl has shown that the corresponding dif- 
fuse scattering is distributed in reciprocal space into a well 
known checkerboard pattern. Figure 1 b) shows cuts in en- 
ergy of this zero magnetic field data set and also of a data set 
with a H ^ 3 T magnetic field applied along the [110] di- 
rection. These cuts simulate constant-Q scans, although they 
integrate the scattering data in the [HHO] direction (shown in 
Fig.l a) and also in the [OOL] direction between L=1.6 and 
1.8. Two features are noteworthy: while there is a shoulder 



to this particular cut of the zero field quasi-elastic scattering 
near E ^ 0.2 meV, there is a continuous and monotonically 
increasing distribution of magnetic spectral weight as the en- 
ergy decreases down to zero energy. The nature of this shoul- 
der near 0.2 meV is subtle for the zero magnetic field data 
set. Indeed, low energy quasi-elastic magnetic scattering is 
obscured by relatively strong nuclear incoherent elastic scat- 
tering which dominates the elastic signal within the 0. 1 me V 
energy resolution of this measurement. However, this nuclear 
incoherent contribution to the elastic scattering can be esti- 
mated and removed by examination of the H — 3T data set 
in Fig. 1 b). Field-induced long range order II20II leads to a 



splitting-off of the quasi-elastic magnetic scattering from the 
resolution-limited nuclear incoherent elastic scattering. Thus 
the strength of the nuclear incoherent contribution can be de- 
termined. Consider the two, otherwise identical cuts shown in 
Fig. 1 b), both at T=0.4 K, and taken al H = (top) and at 
H = 3 T (bottom). One observes scattering at E = with 
an intensity of ~ 9 for H = where the zero energy scat- 
tering has contributions from both nuclear incoherent elastic 
scattering and quasi-elastic magnetic scattering. In the bot- 
tom of panel b), for H = 3 T, an E = intensity of ^ 3.5 
is measured, which has contributions from nuclear incoherent 
scattering alone. Hence the intensity of the elastic magnetic 
scattering at H = is '-.^ 5.5 in the intensity units employed 
in Figs 1 b); at least a factor of two higher than the intensity 
associated with the shoulder near ~ 0.2 meV. Therefore, the 
distribution of magnetic scattering intensity does indeed peak 
at zero energy and extends out to ~ 0.3 meV at low temper- 
atures in H = 0. In other words, there is no obvious mode 
centred at an energy of 1.8 K (0.16 meV) in the H = data 
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FIG. 1: a) High energy-resolution inelastic neutron scattering data 
taken with A=4.8 A neutrons on DCS jioll. b) Energy cuts of the 
data shown in a) as described in the text. Quasi-elastic magnetic 
scattering at T=0.4 K and zero applied magnetic field peaks at zero 
energy. The crystal field excitation in zero field (top right panel) is 
anomalously broad in energy, compared with resolution-limited spin 
waves seen in the bottom right panel within the field-induced ordered 
phase. 

To more definitely assess whether a finite energy excita- 
tion may exist at £' ~ 5 ~ 0.16 meV ^ 0.18 K, very high 
energy-resolution neutron scattering was carried out on the 
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same single crystal of Tb2Ti207, in zero magnetic field, as 
well with several different magnetic field strengths applied 
along the [110] direction. This data, taken with A = 9 A inci- 
dent neutrons on DCS, is shown in Fig. 2. The scattering data 
shown in Fig. 2 integrates the raw inelastic scattering data 
over all Q surveyed, which extends out in reciprocal space to 
the centre of the aforementioned diffuse magnetic scattering 
checkerboard at Q = 002 li, 11|. This very high energy- 
resolution data clearly shows the quasi-elastic magnetic scat- 
tering to increase continuously with decreasing energy down 
to the energy resolution of the measurement, ^ 0.02 meV ^ 
0.2 K, an order of magnitude lower than the singlet-singlet 
gap (5 ~ 1.8 K invoked in Refs. [ilEl. 
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FIG. 2: Very high energy-resolution inelastic neutron scattering data 
employing X = 9 A incident neutrons with DCS, as described in the 
text. The zero field measurement shows the quasielastic magnetic 
distrution of scattering to extend down to at least 0.02 meV. 

Figure 2 also shows that application of a sufficiently strong 
[lIO] magnetic field {H ~ 4 T) moves the quasi-elastic mag- 
netic scattering out of the field of view of the figure, revealing 
a low background which is energy-independent. As a function 
of increasing applied magnetic field, the quasi-elastic scatter- 
ing is strongly depleted at low energies, and an inelastic fea- 
ture is observed for field strengths above H ^ IT. This exci- 
tation moves to higher energies with increasing field strength. 
However, and crucially significant for our argument, we see 
that even with an applied field of H = IT (red curve), the 
broad inelastic peak is at an energy of ^ 0.1 meV ^ 1 K, 
roughly a factor of 2 lower in energy than the gap (5 ~ 1.8 K 
invoked in Refs. |18, 3 in zero field. 

One could speculate that the quasielastic distribution of 
magnetic scattering in zero field out to ~ 0.3 meV is a disper- 
sive singlet-singlet excitation with a band width of twice the 
mean separation between the two states, such that its density- 
of-states fills in th quasi-elastic energy range. However, the 
intensity would then be the largest for the top of this band, 
as the density-of-states would be high where the dispersion is 
flat. This is not observed - the quasi-elastic scattering (Fig. 2) 
decreases monotonically with increasing energy. 

One could also speculate a static and random Jahn-Teller 
distortion in which a broad distribution of singlet-singlet gaps 



is present in the spin liquid state. This could lead to a dis- 
tribution of gaps and induce monotonically-decreasing (as a 
function of energy) quasielastic scattering in zero field as seen 
in Figs. 1 and 2. However, with decreasing temperature, 
these static gaps would get progressively frozen out, result- 
ing in a larger and larger fraction of the system being in a 
non-magnetic state. The l/Ti muon spin relaxation rate is 
large and flat in temperature from ^ 1 K down to 0.05 K Q]- 
This would seem to rule out this scenario since, at 0.05 K, ex- 
cited singlets with an energy gap larger than 0.005 meV would 
be frozen out. Similarly, neutron spin echo (NSE) results in 
Ref. 10] rules out a large 1.8 K gap in the system. NSE with 
sub-/ieV resolution and polarised neutron diffraction |0| both 
show an increase in the magnetic scattering below 300 mK 
suggesting a magnetic ground state or an extremely small gap 
that is thermally active at 50 mK. 

Taken altogether, we conclude there is no compelling ev- 
idence for a well-defined singlet-singlet gap in Tb2Ti207 in 
zero field at low temperatures. Its quasi-elastic, magnetic 
spectrum is not substantially different from that displayed by 
Ho2Ti207 112 in just above its frozen spin ice ground state, al- 
beit with a higher energy scale. 

Residual entropy — Having established that there is no di- 
rect evidence for a split doublet with an energy gap larger than 
0.02 meV ~ 0.2 K in Tb2Ti207, we now address the inter- 
pretation given in Ref. lITsIl of low temperature magnetic en- 
tropy data determined from the heat capacities of Tb2Sn207 
and Tb2Ti207. The magnetic entropy results for Tb2Sn207 
are given in Fig. 2 of Ref. llST . In this data, the nuclear 
and phonon contributions have been subtracted. One sees that 
roughly Rln 4 is lost upon cooling the sample from 20 K down 
to 40 mK. This is consistent with having no extensive residual 
entropy and with magnetic ions having four states lying be- 
neath 20 K. The data for Tb2Ti207 also indicate a similar loss 
of magnetic entropy of Rln 4 below 20 K IT^ . The authors 
of Ref. llSfl proceed to make the following argument: Sup- 
pose that there is a doublet-doublet crystal field scheme with 
a gap A. Then, at high temperatures exceeding the scale of 
the gap, the entropy would be Rln 4 whereas at low tempera- 
tures, T <IC A, the entropy would saturate at Rln 2. There- 
fore, in this scheme, the total entropy variation should be 
Rln 2. Since the entropy variation is observed to be Rln 4 in 
both Tb2Sn207 and Tb2Ti207, Ref. H concludes that the 
doublet-doublet energy level scheme must be incorrect. 

In an attempt to support their argument, the authors present 
a calculation of the entropy for a crystal field scheme consist- 
ing of a pair of low-lying singlets separated by a tuneable en- 
ergy gap of 5 and one excited doublet at energy A > (5 above 
the ground singlet. Upon increasing 5 from zero, the lowest 
temperature entropy drops to zero for T < 5 and exhibits a 



plateau at Rln 2 for 5 < T < A (see Fig. 3 of Ref. UJ]). The 
width of the plateau decreases as T increases and for (5 « 1.8 
K, the entropy variation for this model roughly matches that 
for Tb2Ti207 (the fit is shown in Fig. 4 of Ref. iH). This is 
presented as evidence for a singlet-singlet crystal field scheme 
in Tb2Ti207. A similar conclusion is reached for Tb2Sn207. 
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However, the authors of Ref. lIlSll have not ruled out the 
possibility that the magnetic entropy is lost through the ef- 
fects of interactions and the concomitant build-up of corre- 
lations as the temperature is decreased below 20 K. This is 
particularly evident in Tb2Sn207 which exhibits a phase tran- 
sition at 0.87 K to a long range ordered phase. This implies 
that the residual entropy for T < 0.87 K K should have no 
extensive contribution, as appears to be borne out by the anal- 
ysis of the specific heat capacity data in Ref. llSll . It follows 
that one must account for the details of the transition in order 
to extract information about the single ion level structure from 
entropy data obtained from the magnetic heat capacity while 
no attempt was made to carry out such an analysis in Ref. lIlSll 
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FIG. 3: Panel (a) shows the heat capacity, Cmag (black) and entropy, 
Smag (red), for a one-dimensional spin 3/2 Ising model with J — 
(solid lines) and J/A = 0.05 (dashed lines). The plateau in the 
entropy at low temperatures is at In 2. Panel (b) shows the entropy 
for a one-dimensional spin 3/2 Ising model with J/A from to 0.1. 
The plateau gradually becomes smaller as J/A increases. 



We now turn to the case of Tb2Ti207, which does not ex- 
hibit a transition to long range order down to the lowest ob- 
served temperature of 50 mK Nevertheless, interac- 
tions are significant in this material as can be inferred from the 
Curie- Weiss temperature [iSJ, the paramagnetic diffuse scatter- 
ing QSSlTtl and the dispersion on inelastic neutron scatter- 
ing peaks for example 141 la, 11211 ■ Indeed, the overall scale of 



exchange interactions in Tb2Sn207 and Tb2Ti207 have been 
estimated to be of similar magnitude lll2ll . Although the in- 
teractions in Tb2Ti207 do not manifest themselves as a phase 
transition, they potentially have a large effect on the entropy 
as we now illustrate with a toy model. We consider a spin-3/2 
Ising model on a chain with ferromagnetic nearest neighbor 
exchange J ( J > 0), 



Of importance for our argumentation, this model has (i) a 
ground state doublet, (ii) no phase transition to long range 
order at nonzero temperature and (iii) a paramagnetic entropy 
of Rln4 similarly to Tb2Ti207. This model can be solved 
exactly by splitting the partition function trace into a prod- 
uct of transfer matrices. Consider first the case J = 0. The 
specific heat Cmag shows a single peak (Schottky anomaly) 
corresponding to the loss of entropy, S'mag, as the excited 
S"^ = ±3/2 levels are depopulated upon cooling [solid lines 
in top panel of Fig. 3(a)]. The entropy exhibits two plateaus 



- one at low temperature (Rln2) and another at high tem- 
perature (Rln4). This is the analogue of the non-interacting 
doublet-doublet model considered in Ref. [[isll . 

Switching on the exchange J causes FM correlations to 
build up. There is no transition, but there is a peak in Cmag 
beneath the Schottky anomaly at T / J ^ 0.1. This is re- 
flected in the magnetic entropy, S'mag, dropping below the 
Rln2 value. The Rln2 plateau gradually shrinks as J in- 
creases [Fig. 3(b)]. We conclude that the disappearance of 
an Rln 2 entropy plateau can occur when interactions are in- 
troduced even in the absence of a phase transition and, there- 
fore, cannot be attributed definitively to a splitting of the dou- 
blet ground state without interactions. 

Conclusion — We have reconsidered the scenario advanced 
m Refs. ill [Hi whereby long range order for Tb2Ti207 is 
evaded because the single ion crystal field states are split in 
zero field, ostensibly due to a tetragonal distortion breaking 
the Fd3m symmetry. We have presented high energy reso- 
lution neutron scattering data to re-examine the case for a 
singlet-singlet splitting, finding no evidence for excitations 
that would indicate a splitting greater than 0.2 K. We have 
also argued that measurements determining the residual en- 
tropy cannot be used to draw conclusions about the single ion 
spectrum without considering the build-up of short-range cor- 
relations at low temperatures. In conclusion, the nature of the 
low temperature state of Tb2Ti207 remains a remarkable and 
unsolved problem in the field of frustrated magnetism. One 
may anticipate further progress in light of the constraints that 
we and others are placing on possible scenarios to explain the 
spin liquid behavior in this material. 
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